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Abstract 



We consider the dynamics of bubble growth in the Minimal Standard Model at 
^ . the electroweak phase transition and determine the shape and the velocity of the 
phase boundary, or bubble wall. We show that in the semi-classical approximation 
the friction on the wall arises from the deviation of massive particle populations from 
thermal equilibrium. We treat these with Boltzmann equations in a fluid approxi- 
mation. This approximation is reasonable for the top quarks and the light species 
while it underestimates the friction from the infrared W bosons and Higgs particles. 
We use the two-loop finite temperature effective potential and find a subsonic bubble 
wall for the whole range of Higgs masses < m# < 90GeV. The result is weakly 
dependent on m#: the wall velocity v w falls in the range 0.36 < v w < 0.44, while 
the wall thickness is in the range 29 > LT > 23. The wall is thicker than the phase 
equilibrium value because out of equilibrium particles exert more friction on the back 
than on the base of a moving wall. We also consider the effect of an infrared gauge 
condensate which may exist in the symmetric phase; modelling it simplemindedly, we 
find that the wall may become supersonic, but not ultrarelativistic. 



1 Introduction 

It has gradually become clear that the baryon asymmetry of the universe may have 
been created at a first order electroweak phase transition. The most important ingre- 
dient, baryon number violation, was demonstrated in the Standard Model by t'Hooft 
IjH and was later shown to proceed rapidly at high temperatures [§, §]. The neces- 
sary departure from thermal equilibrium is supplied by the first order electroweak 
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phase transition ||, which has recently been the subject of intense investigation 
10- 0> 0> H> B H9 HID- There have also been promising developments in understanding 



the mechanisms by which baryogenesis may proceed. 

The general picture is this: at high temperatures, such as those prevalent in the 
early universe, thermal effects prevent the breaking of electroweak symmetry by the 
Higgs mechanism. As the universe expands and the temperature falls, the plasma 
supercools in this "symmetric" phase until the probability of bubble nucleation is 
large enough that bubbles of the lower temperature "asymmetric" phase, in which 
a Higgs field condensate breaks electroweak symmetry, form. These bubbles are 
thermodynamically favorable, so they expand, converting the symmetric phase into 
the asymmetric phase. At the bubble surfaces, where the phase conversion occurs, 
the plasma is thrown out of equilibrium by the motion of the phase boundary, or 
bubble wall. Inside the bubbles, baryon number is approximately conserved; outside 
it is rapidly violated. If the departure from equilibrium on the bubble surface biases 
the rate of baryon number violation in a CP violating fashion, this produces a baryon 
asymmetry. 

Specific mechanisms address how the motion of the bubble wall biases the baryon 
number violating processes. The most efficient mechanisms rely on transport. The 
bubble wall separates particles and their antiparticles in a CP violating fashion. In 
the thin wall case this occurs through CP violating quantum mechanical reflection 
|T2| , |13| , |I4| , while in the thick wall case the energy levels of particles and their an- 
tiparticles split in the presence of a CP violating condensate, leading to a different 
perturbation in the population densities of particles and antiparticles |15|, [R| |17| . 
The difference between particle and antiparticle populations is then transported into 
the symmetric phase, where it biases the baryon number violation rate. The models 
based on transport, like most, depend intimately on the details of the bubble wall 
shape and its motion. If the bubble wall is very thin, quantum mechanical reflection 
is the correct language; if it is very thick, scatterings are frequent, and the problem 
should be viewed semi-classically. Correctly describing the situation where the wall is 
thin in comparison to the de Broglie wavelength of some particles, but thick enough 
that particles scatter frequently on the wall, is still an open problem. In the thin wall 
case the main effect comes from the low-momentum particles while in the thick wall 
case the effect from thermal particles dominates. The efficiency of transport strongly 
depends on the velocity and thickness of the wall. For instance, for a slow wall, 
particles diffuse far in front of the wall, and nonlocal baryogenesis is possible. For a 
thin and supersonic wall a typical reflected particle travels about one diffusion length 
ahead of the wall, while for a thick and fast wall no forward transport is possible. In 
the latter case only local baryogenesis occurs, i.e. CP violation and baryon violation 
take place at the same point in space. In any case it is clearly of crucial importance 
to know the shape and velocity of the expanding bubble wall. 

While there has been substantial progress on this problem [|18|, [5], [19], |20|, |21| , much 
remains to be done to compute the wall velocity and shape reliably. In particular, no 
work fully addresses the following issues: (a) infrared boson populations; (b) effects of 
transport; (c) systematic treatment of all relevant scattering and decay diagrams to 
leading order; (d) change in velocity and temperature of the plasma due to latent heat 
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release; (e) dynamical determination of the wall shape (thickness). In a recent letter 
we developed a set of techniques and approximations which address (b), (c), and 
(e). In this paper we present the details of that work and extend it to include (d). 
The problem of infrared boson populations is of theoretical interest, since it is related 
to the infrared properties of gauge theories at finite temperature, which are not well 
understood. We leave the proper treatment of this problem to a future work. The 
fluid approximation we use probably underestimates the friction from infrared bosons, 
so our result could be interpreted as an upper bound for the wall velocity. Also note 
that we work within the Minimal Standard Model, even though extensions with more 
CP violation are generally considered more viable candidates for baryogenesis. We 
do this for simplicity, developing the techniques to control the problem before we 
attempt to apply them to more complicated and interesting models. 

Now let us briefly outline the paper. In section II, we derive a semi-classical 
equation of motion for the Higgs condensate in the presence of out of equilibrium 
particle populations. The semi-classical approximation is accurate provided the scale 
on which the Higgs field condenses, given by the wall thickness, is large in comparison 
to the inverse momentum scale of particles which significantly contribute to friction. 
Solving this equation requires knowledge of the effective potential, the temperature at 
the wall, and the departure of particle populations from equilibrium in the presence 
of the wall. In section III we review the state of the art in the effective potential. In 
section IV, we study the bubble nucleation using Langer's formalism and the potential 
of section III. This formalism was applied to the field theory by Coleman and Callan 
23] and to the electroweak phase transition in |24], |, [19], and permits the 
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calculation of the temperature at which the bubbles nucleate. 

The liberation of latent heat produces a jump in temperature and velocity across 
the wall, and changes the temperature behind the wall with respect to the temperature 
far from the bubble. We compute this effect using stress-energy conservation, as 



done for example in p4| , |19| , p6| . Including this effect could be important since it 
may significantly increase the friction on the wall, especially when the wall velocity 
approaches the speed of sound. This is an important improvement on our work in 
pifl , where for simplicity we ignored this effect. 

Next, we address the departure of particle populations from equilibrium in the 
presence of the wall. In section V we solve exactly a particularly simple case of free 
particle scatterings off the wall (assuming efficient diffusion so that no piling up occurs 
in front of the wall). Unfortunately this is a reasonable approximation only for an 
unrealistically thin wall. However, it does help us to determine what aspects of the 
departure from equilibrium are most important. For fermions we find that most of 
the friction arises from thermal energy particles, while for bosons most arises from 
infrared particles. The friction depends on a particle's mass as m 4 , so top quarks 
appear to be the largest contributors to the friction, even though they are scattered 
fairly efficiently This justifies the development of an approximation for the particle 
populations which models thermal particles reasonably accurately. This is the fluid 
approximation, which is developed in some detail in section VI, with a systematic 
study of the leading order tree-level processes that restore equilibrium postponed 
to Appendix A. Section VII studies aspects of transport in the fluid equations by 
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finding the response to a ^-function source. It is shown that, in contrast to a slow 
wall, no forward transport is possible for a supersonic wall. The accuracy of and 
possible improvements to the fluid approximation within the context of a momentum 
expansion are discussed in Appendix B. 

Finally, we combine the fluid approximation with the equation of motion for the 
Higgs field and solve for the bubble velocity and profile. In spite of the nonlinear 
character in the equation of motion, this system of equations can be solved analytically 
using the Fourier transform and a two parameter Ansatz for the wall shape, as we 
show in section VIII. In section IX we describe a numerical technique which allows 
a general wall profile. The results are presented in section X. The numerical work 
shows that the Ansatz models the wall velocity with good accuracy (typically within 
5%) although it is not as good at modelling the wall profile, as illustrated in Figure 
|3|. We conclude that the Ansatz technique is reasonably accurate for subsonic wall 
velocities. However, for very supersonic velocities, which may occur if infrared physics 
generates gauge condensates in the symmetric phase, both the Ansatz, and the fluid 
approximation itself, break down. This is discussed in Appendix C, where it is shown 
that the bubble wall cannot propagate as an ultrarelativistic detonation. 

For the impatient reader we suggest a fast track to reading the paper: section II, 
the second half of section IV, section VI, and sections VIII - X. 

2 Equation of motion 

We are interested in the dynamics of an infrared Higgs condensate, which we will 
treat as a classical background field. From the terms in the Electroweak Lagrangian 
containing $, 

where the sum is over all massive fermions and y is the Yukawa coupling, we derive 
the operator equation of motion 

2 

□$ + igA^d^ + y (d ■ i)$ - ji 2 $ - //$ + 2A(<f> t $)$ + V$r&l 

We have suppressed group indices. Each term containing A actually appears once for 
weak isospin and once for hypercharge. 

This operator should annihilate the physical thermal state. We will shift $ by 
a classical part, $ = & c i + <5$, choosing 5$ such that (5$) = 0. We then evaluate 
the trace of the operator equation of motion over the (out of equilibrium) thermal 
density matrix describing the propagating bubble wall. We will assume that there 
are no charge conjugation violating gauge condensates, so that (A^) = 0. We then 
find 

□ $ d - + 2A($t / $ d )$ c ,+ 
2A (2(6&6*)* d + (5$ 2 )$y - ^-(A 2 ) + £ y@ R if> L ) = 
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We assume that $ c z = [0 0/ y/2 } T (which is the same as neglecting charge conjugation 
violating condensates). Next, we evaluate the thermally averaged operators using 
WKB wave functions. This makes sense because the background field varies on a 
scale much longer than T" 1 , which characterizes the reciprocal momenta of particles 
in the plasma. In this approximation, 



where the sum is over group indices and spins and / is the phase space population 
density, and 

l „ ^ r d 3 k 



where the sum is over group and Dirac indices. 
The equation of motion becomes 



□, + n , ) + E !g/_**L_ /(M) = o (1) 

where V is the renormalized vacuum potential, and the sum includes all massive 
physical degrees of freedom, including the high frequency parts of the Higgs field. 
Note the condensed notation: m 2 = y 2 <j) 2 /2 for quarks and leptons, <? 2 2 /4 for the 
gauge fields, 3A0 2 — /x+ (thermal) for the Higgs bosons, and A0 2 — /i+ (thermal) for 
the pseudo-Goldstone modes. This equation has also been derived by diagrammatic 



techniques in 19 



This equation has a nice physical interpretation. Multiplying it by <9 M 0, we find 

f d k 

= D0^0 + d»V + J ——f(k)d»m 2 (2) 

The first term can be recast as d^d^(j)d v (j) = ((9^0(9^0) — d u {d^ l (j)d ll (j)/2) , so that 
when taken together the first two terms are the divergence of the stress-energy of 
the Higgs field T^(0) = d^(j)d v (p - £((/>), where £(0) = ^0^0/2 - V(<j>) is the 
Higgs lagrangian. The last term represents exchange of stress-energy with particles; 
it looks like J d 3 kf(k)(-F^)/(2ir) 3 . Here F» = -d"E = (-d t E, -VE) is the 4-force 
a particle feels in the presence of the wall; the wall feels an equal and opposite force. 
Thus we can recast (0) as follows 

/d k 
— /(fc)F" = 0, FP = -d»E (3) 

The equation of motion for the Higgs field in the WKB approximation splits the total 
(conserved) stress-energy tensor T^ u into Higgs and fluid parts (sub-systems) in a 
natural way, and provides the prescription for how each is violated. 

As an application of (|3|), J*(Eq.|]) <p'dz is the total pressure on the wall. If this 
integral does not vanish, then the wall accumulates momentum and accelerates; hence 
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this integral must vanish for a steady state wall. This fact forms the basis of the 
analysis of [0, |, [BJ. 

We write the population density / as the equilibrium population f plus a devi- 
ation, / = /o + 5f. The vacuum contribution V'(<f>) and the contribution from fo 
combine to give the finite temperature effective potential V^(<f>). Thus we have 

^ + vW) + £^/p|^/fe*Ho (4, 

If the system were in phase equilibrium we would have expected □</> + Vj.; we see that 
the frictive force (dissipation) arises due to the departure from thermal equilibrium 
5f (fluctuation). Our goal will be to solve this equation for a stationary wall, i.e. 
after the wall has reached a steady shape and velocity. To do so we need to know the 
effective potential Vr and the temperature T, and we need a way to calculate 5f. 



3 Effective Potential 

The high temperature expansion of the one loop effective potential, ignoring scalar 
loops, is l| 



Here 



V(0, T) = D(T 2 - T 2 )0 2 - ET(j) 3 + . (5) 

D = T^(2m 2 w + m| + 2m 2 ) ~ 0.169 , 

E = ^} 2m w + m%) ~ 10- 2 , 
To=^«-8Bv 2 ) , 

• (6) 

where lna B = 21n47r — 27 ~ 3.91, lna^ = 2 In ir — 27 ~ 1.14, and 

5 = 7^rr( 2 ^ + <-4m 4 ) 



647T 2 f 



The value of in the broken phase at equilibrium is easily evaluated to be 2ET / At- 
If we take At parametrically ~ g 2 , which is natural from the renormalization structure 
of the standard model, we find <fi ~ gT, which is small enough that the perturba- 



tion series may not be well behaved |27[| . We should therefore include the two loop 
contribution, which has recently been computed by and . The two loop expres- 
sion contains terms ~ g 4 <p 2 T 2 , with a coefficient including constants and logs of the 
temperature over the renormalization scale, which slightly corrects D. We will drop 
this term as it has no influence on the behavior of the phase transition. Two more 
important corrections do occur, however. One is that the longitudinal components 
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of the gauge bosons acquire large plasma masses ]28|, |5| and do not contribute to 
the strength of the transition. Following ]5| we drop them from the E term in the 
effective potential, giving 

E = -^(im 3 w + 2m 3 z + (3 + 3 L5 )A^ 5 ) 

where the At dependent term comes from including ring improved scalar loops. The 
other correction is that qualitatively new terms ~ g A <f) 2 T 2 ln(m/T) appear. Here m 
stands for sums of masses of particles. In the approximation that the Higgs masses are 
small in the symmetric phase (which should be good near the transition temperature) 
we may treat m ~ g<p; the log of g is absorbed into D, and the new term becomes 
-C(f) 2 T 2 \n((j)/T). The value of the coefficient is 

C ~ ^— 2 {lA2g A w + 4.8(7= A T - 6A T ) 

where we have again left out contributions from transverse gauge bosons. Including 
this term, the value of <fi in the broken phase at equilibrium becomes <fi/T = Ej At + 
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\t) 2 + 2C/Xt, which indicates that the new term strengthens the transition 
and contributes at the same parametric order as the one loop term E, although for 
the parameters we will be interested in its contribution is smaller. 

Because the two loop result is of the same parametric order as the one loop re- 
sult we might worry that perturbation theory cannot establish the strength of the 
transition reliably. In fact we expect that the perturbative computation of the ef- 
fective potential should break down in the symmetric phase. This means that the 
value of Vt(0) may be shifted somewhat from zero. Shaposhnikov has suggested that 
such a shift may arise due to the formation of gauge condensates in the symmetric 
phase f8|. He advocates adding a term —Apg^jYl * Pit(0) to the effective potential 
to account for such an effect. Here the function Pit describes the </> dependence of 
the strength of such a condensate, and is smooth at <fi = but falls as exp(— <f)/g 2 ) 
at larger 0. We will include such a contribution to parametrize our ignorance of the 
effective potential in the symmetric phase. We choose Pit=sech(0/ayt/) because it is 
simple, but our results do not depend strongly on the functional form or the exact 
rate of exponential decay, as long as that decay is rapid. We treat the value of At 
as an unknown parameter; there is some lattice evidence for its value but it is still 
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Our final form for the effective potential is then 

V T {4>) = D(T 2 - T o 2 )0 2 - CT 2 <f In A - ET<f + ^ - sech( A) (7) 

4 Nucleation and Hydrodynamics 



We must next determine what T is at the bubble wall. Following 0] we do this 



in two steps; first we calculate the temperature in the universe when most bubble 
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nucleation events occur, and then we calculate the temperature correction due to the 
release of latent heat by the bubble as it expands. 

The process of bubble nucleation at finite temperature was worked out by Linde 
and applied to the electroweak phase transition in [18], ^], [19], |25|. The 



basic idea is that small bubbles are thermodynamically unfavorable because of positive 
surface free energies; to become large enough not to re-contract a bubble must pass 
over a free energy barrier. The lowest route over the barrier goes through a saddlepoint 
configuration of the effective action called the critical bubble. It is a spherically 
symmetric solution to the classical equilibrium equation of motion 

V 2 0(r) = V+(<t>) <f>'(r = 0) = , 0(oo) = (8) 

The free energy of such a solution is 

Scrit = / rf 3 xQ(V0) 2 + K r (0)) (9) 

There are analytic formulas for S CT i t only in special cases; in general it must be 
determined numerically. 

The nucleation rate per unit volume per unit time is 

I = / T 4 exp(-^) (10) 



Carrington and Kapusta have recently computed the prefactor ||25|| . Numerically eval- 
uating the expressions in their paper we find, roughly, ln(/ ) — —14. The exact value 
depends on the parameters of the effective potential and (weakly) on temperature, 
but as we will see we only need a rough estimate. 

We can now compute the action of the critical bubble at the time when most of 



the universe changes phase by following the technique of |3(J • If bubbles expand at 
a velocity v w which depends only weakly on temperature, then the fraction of space 
remaining in the symmetric phase at time t is roughly 

exp (- fvl(t - t')%e- s ^ T dt') (11) 



where the exponential accounts for bubble overlap ||31|| . Expanding the bubble action 
about the point when the above integral is 1, S = S(T nnc ) + (t — t')dS/dt and using 
d/dt = (dT/dt)d/dT and (l/T)dT/dt = H the Hubble constant, we find that most of 
space has changed phase when 

5(T nuc )/T _ 87r ^ J 

(HTdS/dTY 

A reasonable estimate for TdS/dT is the thin wall value 2ST/ (T eq — T) ~ 10 4 , which 
roughly agrees with the value we find by numerically evaluating S/T at closely spaced 
values of T. The value of H follows from the Friedmann equation 

3 30 
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where the second factor is the energy density of the plasma and g>* = 106.75(1 +0(a s )) 
in the symmetric phase (the 0(a s ) arises from thermal masses and other interaction 
effects). The value of v w is to be determined, but since a subsonic bubble is proceeded 
by a hydrodynamic front travelling at the speed of sound which raises the temperature 
of the fluid and prevents nucleation, it is sensible to choose the speed of sound v s = 
1/V3. 

Putting all the expressions together and using T ~ lOOGeV, we find S/T ~ 97. 



In Ref. |MJ a slightly higher value was found because of the value for I we use. We 
may now determine the temperature at which most nucleation occurs by solving (|) 
for various values of T until we find one which gives S CI n = 97. 

This gives the temperature around the bubble at the time that it nucleates, but we 
are really interested in the temperature at the bubble wall while it is expanding. This 
will be elevated from T nuc because of the liberation of latent heat as the symmetric 
phase is converted into the asymmetric phase. 

In H], |1| it is argued that the elevation of temperature at the bubble wall is 



not important in determining its velocity. This is true when the change in temper- 
ature, which is ~ l/(dp/dT) (where I is the latent heat), is small compared to the 
supercooling, T eq — T nuc . [19|] show explicitly that for the effective potential they 
consider, (T cq — T nnc )(dp/dT) is about 5 times /. However, that analysis is based on a 
small assumed value of m t ; the latent heat depends onDoc m 2 , and the supercooling 
T eq — T Q oc 1/D; so for m t ~ 174GeV, the temperature increase probably is important. 

To compute the temperature at the wall we need to solve the temperature profile 
around an expanding bubble wall. The equation of state in the symmetric phase is 
p s = (7r 2 (7*/30)T 4 , p s = p s /3, and in the broken phase it is p a = p s + /(T), p a = 
Vs ~ Vt(0(T)). Here I is the latent heat and is given by 

I = W«) - T d ^ L = ~D<P 2 (T* + T 2 ) + - 5)0 4 ~ -2D0 2 T 2 (12) 
dl 4 

Note that l/p ~ 0.01 << 1, and similarly we find (T eq — T)/T ~ 0.01 << 1, so it is 

reasonable to expand to lowest order in these quantities. However, / turns out to be 

quite a strong function of temperature, 

d(Z/T 4 ) 6Td(6/T) , N 

dT T dT K J 

We can determine Td(<p/T)/dT in the broken phase as follows; the value of <ft is 
determined by V^((p) = 0, and since nonperturbative effects are tiny in the broken 
phase we can use the two loop expression, which gives 

= 2£>(T 2 - T 2 )0 - C0 2 T 2 (2 In j- + 1) - 3£0 2 T + A T 3 

= 2D(l-S)_2Cln^-C-3^+A T ^ (14) 

Since this equals zero, its total derivative with respect to T, d/dT + d(<p/T)/dT 
d/d((j)/T), should vanish. This gives 

d{cj>/T) -ADT^/T 



dT ~ 2X T (p 2 - 3E(pT - 2CT 2 



(15) 
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and hence 

d(l/T 4 ) 4D<P -ADT^/T 

dT T 2X T (f) 2 — 3E(j)T — 2CT 2 1 ' 

For typical values of At = 0.03,0 = 0.8, we find Td(l/T 4 )/dT ~ 18, which is about 
half as large as p/T 4 . This rather surprising result occurs because <fi is strongly 
temperature dependent in the broken phase near the transition, and because the top 
quark is so massive; around the transition temperature the number of top quarks is 
rapidly freezing out. Because p turns out to have this extra temperature dependence, 
the speed of sound in the broken phase, dp/ dp = (dp / dT) j (dp / dT) , is around 15% 
smaller than in the symmetric phase. Thus it is reasonable to expand to linear order 
in l/p and (T — T eq )/T, and hence in ST and v , but not to neglect dl/dT. 

Now let us explore the fluid temperature and velocity in the vicinity of an isolated 



spherical bubble. From Eq. ( |TT|) and the argument which follows, we find that the 
typical bubble grows to a radius of ~ 1/ (HTdS/dT) ~ 10 12 /T, while the distance it 
takes for a bubble to accelerate to a steady velocity is around the ratio of the surface 
tension to the free energy density difference of the two phases, which is of order 50/T. 
On the scale of the bubble radius, then, we can treat the plasma as a perfect fluid 
and expect the temperature and velocity to be functions of r/t = x only, where time 
t is measured from the time of nucleation. Conservation of the stress energy tensor 
gives the equations 

= d^ v = d^((p+p)u^u u -pg^) = u u u-d(p+p) + (p+p)(u u d-u+u-du u )-d u p (17) 

which to linear order in v and ST (which are linear in l/p) give 

, / 2 ^ dv\ dp dT g dp dT ^ . dv 

\x dx J dT dx dT dx dx 

The 2v/x arises from taking a divergence in spherical coordinates and is the only 
difference between the spherical geometry and the planar geometry studied by p4| , [19]. 

Because the equation of state for the symmetric phase is simple we can solve the 
above equations for the case that the wall velocity is less than the speed of sound in 
the broken phase (in which case T is constant and v is zero inside the bubble). The 
equations become 

dv d\nT ( x 2 \ dv n 

x— = — — 2v + x[l- — — = (19) 



which give 



dx dx 



^ 1 1 + -) (i - -) to ^ K i^- = 2A / s.) (21)) 

X / \ X / 1 niir J- nur. V X 



Here A is an undetermined parameter and v s = l/y/3. It is interesting to note that, 
in this approximation (linear order in l/p), the temperature and velocity are actually 
continuous at x = v s , although the first derivative is discontinuous. This holds to all 
perturbative orders in l/p. The shock front is exponentially weak, as found in [|(| . 
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Since the wall does not accumulate energy or momentum, Eq. ( |18]) also provides 
us with boundary conditions at the bubble wall, by replacing the derivatives with 
differences. Because the wall is thin compared to the bubble radius the term 2v/x is 
negligible and we get 

(Ps + Ps)v - V w p s = -V w p a (T inside ) - V w (p s + p s )v +p s = p a (T inside ) (21) 

which, at fixed v w , allows us to solve for the temperature in front of and behind the 
wall. Defining y = v w /v s , we find 

T s -T nuc = y 2 l(T a ) T a - T nuc y 2 l(T a ) 

T nuc 2(l + y)2(l-y) T nuc 4(1 + y) 2 1 j 

Note the dependence of I on T a the temperature inside the bubble, which makes these 
equations transcendental. 

When the wall velocity exceeds v s , the only causal solution to ( ]T9| ) outside the 
bubble is v = 0, T = T nuc . We may find the temperature just inside the bubble wall 
directly from the boundary conditions 

(pa + Pa)v - v w p a = -v w p s (T nuc ) - v w (p a + p a )v + p a = p s (T nuc ) (23) 

which give 

T a -T nuc _ y 2 l(T a ) 



T nuc 4(y + l)(y-l) (24) 

In both cases the temperature is elevated on one or both sides of the bubble. The 
elevation has a l/(y — 1) behavior, which diverges near v s . This divergence means that 
our expansion in small v and T — T nuc must break down very near v s ; but there is still 
a very large elevation in temperature when v ~ v s . As the symmetric phase becomes 
more favorable at higher temperatures this prevents the wall from propagating close 
to the speed of sound; it must either be a deflagration or a strong detonation. 

We have not treated the case where v w lies between v s and the speed of sound 
in the broken phase; this case is much more complicated because the fluid velocity 
is nonzero on both sides of the wall and we need three conditions to determine all 
unknowns. Fortunately, in this case the temperature would be elevated enough that 
the wall probably cannot propagate at this speed; so we will not treat this case. We 
will also not treat the effect of latent heat released by other bubbles, although as 
the phase transition progresses and bubbles begin to collide, one bubble's wall will 
frequently be in the hydrodynamic wave initiated by another bubble. These effects 
make the later phases of the transition quite complicated, and the wall velocity we 
compute here should be an upper bound on the velocity in this epoch p6[ . This effect 
should be important especially when the wall velocity is much smaller than the speed 
of sound. In this case bubbles would expand most of the time in the background 
of a sound wave of another bubble. This effect may significantly enhance baryon 
production clS duT erued in IE 
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5 Free Particle Approximation 



Let us proceed to computing the wall velocity. First we will discuss a simple limit, 
the case where mean free paths are enormously longer than the wall. This limit has 
been treated thoroughly by [gj, |19| and we will only discuss it briefly to illuminate a 
few points about which particles provide most of the friction. 

In the WKB (semi-classical) approximation particles should obey a Boltzmann 
equation, 

d t f + S-d 3 f + f-d ff f = -C[f\ (25) 

Here p = —zdE/dz = (l/2E)dm 2 /dx, dtf = +v w d z f (where we choose the broken 
phase at z — +00 so the wall moves to the left) and C[f] is the collision integral. 
The free particle approximation consists of dropping C[f}. Since collisions are what 
drive / to equilibrium, we must insert equilibrium by hand as a boundary condition 
for incoming particles. We will use this approximation (though we should note that 
some of the departure from equilibrium may reflect back onto the wall [21 1). In this 
case the Boltzmann equations can be solved exactly, 



/ 1 = exp I ' - ^-f I yr \ ' — ' + m z I \ ±1, p z > - Jm'^-m 



E + v w p z ~/v w ( (p z + v w E) 



2 



(l-v 2 w )T T V l-v. 

, E + v w p z jv w f(p z + v w E) 2 2 ., , , 

6XP ' {l-v*)T ~T~ { l-vl m ~ m ° ' ' ' Pz ~V ~~ ( ^ 

Here m is the mass in the broken phase and ± is + for fermions and — for bosons. 
Determining J d 3 k/((2n) 3 2E)f , the term in Eq. (|J), for this expression appears 
hopeless, but we can expand to lowest order in v w , which gives 




/ 



d 3 p 

Jo + 



(2tt) 3 2£ 



d 3 p e E/T J ^Jpl + m 2 - p z , p z > -v' mr, - - m- 



{2tt) 3 2E {e E / T ± l) 2 I -Jpl + m 2 _ m 2 _ pz ? pz< _/ m 2 _ m 2 



(27) 



This term is independent of the momentum perpendicular to the wall, p±; so we can 
perform the integral over p±. 



d 2 p ± e E / T _ fo^P 2 z + m 2 



(2n) 2 2E {e E / T ± l) 2 4tt 

Integrating this expression times mm' over the wall yields the velocity dependent part 
of the pressure; the result agrees with the integral found by 0, who have evaluated 
it to order m 4 . 

Let us examine the integral over p± more carefully. For bosons, the integral is 

P±dp± 1 2 E 

csch — — 

AttE 4 2T 
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At small E, csch (E/2T) ~ 4T 2 /E 2 and the integral has a linear small momentum 
divergence, cut off by Jp 2 + m 2 ; the dominant contribution is infrared. The resulting 



integral over p z is also linearly IR divergent, as both / and ^Jp 2 + m 2 — p z behave as 
l/E above itlq. Thus, the friction from W bosons arises primarily from very infrared 
(E ~ m) particles. This is because the Bose distribution function has a pole at E — 0, 
and because a particle's contribution to the Higgs equation of motion goes as l/E. 

Fermions have much more mild infrared behavior; the integral over transverse 
momenta is 

p±dp± 1 2 E 

sech — 

AttE 4 2T 

and at small E the function sech 2 (E/2T) behaves as 1; the integral is well behaved in 



the infrared, even for ^p 2 + m 2 = 0, and the dominant contribution is from particles 
with px ~ (1 — 2)T. However, the integral over p z is small p z divergent, though only 
logarithmically; so the dominant contribution from top quarks comes from particles 
with small transverse momenta, but thermal energies, i.e. from particles at glancing 
incidence. This is a combined effect of the large phase space in such particles and 
their long time on the wall. 

Now let us examine the analysis of Dine et. al. |J. They propose to treat the 
particles in a relaxation time approximation. However, they compute the particle 
population at a point as having come off of a sheet a fixed z away from the point of 



interest. A particle's lifetime then actually goes as rp/p z = rJl + (p±/p z ) 2 i so the 
relaxation time is longer for particles travelling at glancing incidence. Their result 
then contains a logarithmic enhancement, which can be traced to these glancing 
incidence particles. 

We will now show that no such log occurs in a relaxation time approximation 
when the relaxation time is short compared to the thickness of the wall. Writing 
/ = fo + 5f and working to lowest order in Sf, and treating the wall as stationary so 
Sf = 8f(z + v w t), the Boltzmann equations become 

[Vw + E }6} + C[n ~ Vw 2E (exp(£/T) ± I) 2 (28) 

The relaxation time approximation is C[f] = 5f/r, with r independent of momentum. 
If r << L (with L the wall thickness) then to lowest order in t/L the derivative term 
may be dropped, giving 

6 f- TV ^l e ^ E ^ (29) 

The friction on the wall from one species in the relaxation time approximation is 
then 

dm 2 r d 3 k (m 2 )' exp(£/T) 

> Z(P ^J (2Kf2E TVw 2E (exp(£/T)±l) 2 (3 ° } 
which, at lowest order in an expansion in m, integrates to 

TV,, 



4tt 2 



/ m 2 (m') 2 dz , (fermions) — — / m{m!) 2 dz , (bosons) (31) 

J 8lT J 
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which for a tanh wall shape give ru^mj/ '40it 2 L and rv m m\j '48ttL respectively. The 
top quark contribution is isotropic, dominated by thermal particles, and has no log. 
The log will be recovered when r >> L and the derivative terms in Eq. fl28|) become 
important. 

We do not expect this approximation to be very good; in particular we should 
include the derivative terms and model C[f] more accurately. This is the goal of the 
remainder of the paper. 



6 Fluid Approximation 

In this section we will present a method for calculating the deviation from equilibrium 
population density 5f in (ffl) in the presence of a moving wall. Our starting point is 
the Boltzmann equation, 

dtf = d t f + zd z f + p z d Pz f = -C[f) (32) 

where C[f] represents the scattering integral, E = (p 2 + m 2 ) 1//2 is the particle energy, 
and 

d Vz . d {m 2 ) 1 

Z = Vz= dp- E= E> Vz = -Yz E = -^E (33) 

are the velocity and the force on the particle, respectively. The Boltzmann equation 
is the semi-classical approximation to the quantum Louville equation. To be valid the 
background field must vary slowly enough that particles satisfy the WKB condition: 

1 

P » -r~ 



L 



w 



Since we expect rather thick wall, L w » 1/T, this relation is satisfied in abundance 
for all but the most infrared particles. We will assume that very infrared W bosons 
are scattered up to energies O(gT) quite readily, so W bosons with E ~ gT or 
larger are responsible for more friction than very infrared W bosons. We will use this 
assumption again and again in what follows. It may be wrongf] , in which case we 
need to extend our analysis with a special treatment of very infrared bosons. This 
extension must model the dynamics of the strongly coupled infrared sector of the 
thermal field theory, which we are not able to do at this time. We will also neglect 
any friction from the condensate responsible for Ap (if such a condensate exists) for 
similar reasons. It is our belief that this problem and the problem of determining the 
friction arising from thermal and O(gT) particles may be treated independently, so 
our analysis may be simply extended once the infrared sector is understood. 

The Boltzmann treatment also requires that scatterings with the plasma are not 
too frequent so that particles can be to a good approximation considered on-shell for 
all times. To quantify this we write the Heisenberg energy-time relation AE ~ 7 



3 The analysis of which neglects transport but attempts to treat the energy dependence of 6f 
more carefully, apparently supports our assumption, but their analysis neglects thermal corrections 
to the infrared boson propagator. 
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where 7 is the appropriate damping rate which is of order g 2 T . Since we can only 
treat particles with E,p > gT this uncertainty can be ignored. 

The Boltzmann equations are nonlinear partial integro-differential relations and 
as such are analytically intractable. To evade this difficulty we will model each 
distribution function with a several parameter Ansatz called the fluid approxima- 
tion. The Boltzmann equations will then yield a set of ordinary differential equa- 
tions for the parameters which will be tractable analytically. The Ansatz is f^ 1 = 
exp(E/T — E5T/T 2 — fi/T — p z v/T) ± 1, the form for a perfect fluid. We should use 
different parameters for each species in the plasma, but we will make an additional 
approximation that all light degrees of freedom (that is, all but top quarks, W bosons, 
and Higgs bosons) are in equilibrium at a common (space dependent) temperature 
T+5Tb g and velocity v bg . We will also give tops and anti-tops of both helicity the same 
distribution function. In the Minimal Standard Model this is reasonable as there is 
almost no CP violation, and the difference in transport properties arises only at the 
subleading level of weak scatterings. In two doublet models we might need to be more 
careful. Similarly, we will treat W bosons and Z bosons as a single species (henceforth 
U W bosons" ) whose mass squared is a weighted average of the mass squared for the 
W and the Z. This ensures that the effective potential parameter D is correct. 

For the heavy degrees of freedom we take the distribution function to be 



f = fo(E + 5) = e(E+s)/T±v 8 = - 



E 

H + Hbg + —(5T + 8T bg ) + p z (v + v bg ) 



(34) 

We track three perturbations with respect to the background; chemical potential /x, 
temperature 5T and velocity v. This Ansatz is a truncation of an expansion in powers 
of momentum; it gives a reasonable description for the thermal particles' populations 
when the wall background varies slowly on the scale of the diffusion length. For top 
quarks this should be sufficient as the diffusion length is short and the influence of 
infrared particles is phase space suppressed. We discuss its validity at some length 
in Appendix B: suffice it to say that it appears quite reasonable for top quarks, and 
quite naive for W bosons. 

We now outline the main steps in the derivation of the fluid equations. First we 
expand d t f to linear order in perturbations 

TdJ = \E + d t S\fL JJ= eME/T) 



(exp(E/T) ± l) 2 

d t 5 = -(^d t + ^d z yii + ^ g + Pz (v + v bg )}-^(d t + ^d z y8T + ST bg ){35) 

This equation is written in the fluid frame so that a particle's energy is not conserved 
in the presence of a moving wall. The term which perturbs the population densities 
away from equilibrium is 

Wo = ^/o (36) 

Note that this term is proportional to m? . This is the reason it is sensible to expand 
to linear order in perturbations. First note that m « T. For thermal particles 
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with E ~ T, — /q ~ /o and the source for Sf/f is 0(m 2 / ET) which is small. For 
infrared fermions, — f' Q ~ /o/2 and the source is 0(m 2 / ET) which is also small. For 
infrared bosons the situation is a little less rosy; at small energy — f' Q ~ foT/E and the 
perturbation is 0(m 2 /E 2 ). Again, we will assume that these particles are scattered 
efficiently and postpone a more careful treatment. 
The Boltzmann equations have now become 

(~fo) (f [dz(» + in*) + ^d z (5T + 5T bg ) + p z d z (v + v bg )\ + d t (fi + fi bg ) 

+^d t (5T + 5T bg ) + Pz d t (v + v bg )) + Ojjl, 5T, v) = (-/o)^f^ (37) 

The collision term depends on the deviations of all particle species from the common 
background temperature and velocity. We discuss it in some detail in Appendix A. 

The three parameters are determined uniquely by taking three integrals of the 
Boltzmann equation; since our perturbations are Lagrange multipliers for particle 
number, energy, and momentum, the appropriate choice is / d 3 p/ (2n) 3 , J Ed 3 p/ (2n) 3 , 
and / p z d 3 p/(2ir) 3 . The resulting equations are 

c 2 d t (» + » bg ) + c 3 d t (5T + ST bg ) + C -^d z (v + v bg ) + / j0f- 2 C[f] = ^d t m 2 

c 3 d t (pt + pt bg ) + c,d t (5T + 5T bg ) + °^-d z {v + v bg ) + / j0j£*C\f\ = ^d t m 2 
C3 , . C4 , ^ , C4T r p z d 3 p 



jd^ + ^ + ^dtiST + ST^ + ^dtiv + ^ + J-^^Clf] = (38) 
The constants q are defined by 

For fermions they are ci/ = ln2/27r 2 , Q/ = (1 — 2 1 ~ t )i\( i /27r 2 , % > 1. For bosons they 
are c\ b = ln(2/m)/27r 2 , Ci b = i\Q/2ir 2 . Here and throughout Q is the Riemann zeta 
function evaluated at %. We have worked to lowest order in m/T here, but only c\ b 
and c 2b possess 0{m/T) corrections, and these are small. 

For a stationary wall all quantities are functions of x = z + v w t and the stationary 
fluid equations are obtained by a simple substitution 

dA^v w 5[, dA ^ 8'i (39) 

The collision integral / C[f] depends on /i, 5T, and v of all heavy species, which 
couples the fluid equations. However, for the top quarks the collisions which con- 
tribute to order a 2 arise only from collisions with light quarks and gluons, and for W 
bosons the contributions at order a s a w and a 2 are primarily with light quarks (the 
top quark making up only 1/6 of the population of quarks). Hence it is reasonable to 
ignore direct collisions between top quarks and W bosons. If we included them then 
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the problem would become more complex, as these interactions mix top and bottom 
quarks, and distinguish between top helicities. To treat them properly we would need 
to introduce separate fluid equations for left handed bottom quarks and for each top 
quark handedness. Because the effect is ~ 10%, we will neglect it and treat all W 
boson collisions with quarks as being with light quarks. 

The Higgs boson interacts predominantly with the top quark via its Yukawa cou- 
pling, but since this interaction is quite efficient (0(a t a s In l/a s )) and the Higgs has 
only one degree of freedom (and is lighter than the W boson if the phase transition 
is to be strong enough for baryon number to be conserved after its completion) we 
will ignore Higgs particles altogether. Hence the collision terms in the fluid equations 
may be treated as arising entirely through interactions with the background of light 
particles. 

We have computed the collision integrals appearing in Q3"5|) , including all diagrams 
which contribute to order a 2 s In 1/ a s for top quarks and to order a 2 In 1/ a w for W 
bosons, in Appendix A. All diagrams are evaluated in the leading-log approximation. 
This means that only t-channel processes (which are logarithmically divergent in the 
limit of the zero exchange particle mass) are calculated. The result for top quarks is 



d 3 p 



(2tt) 3 T 2 
d 3 p 



C[f] = fiT^f + 6TT Tlf r Ml/ = .00899T T Tlf = .01752T 



(2vr) 3 T 3 



EC[f] = fiT/af + 5TT T2f T^f = .01752T T T2f = .07172T 



d 3 p 



(2tt) 3 T 3 



Pz C[f] = vTT v f rV = .03765T (40) 



We use the subscript / (fermion) for top quarks and b (boson) for W bosons, but we 
suppress them when no confusion will occur. 
For W bosons we get 



d 3 p 



(2tt) 3 T 2 
d 3 p 



C[f] = + 6TT Tlb = .00521T 



(2tt) 3 T 3 



EC[f] = fxT^t + 5TT T2b r M2b = .01012T 



T Tlfe = .01012T 
T T2b = .03686T 



d 3 p 



Pz C[f}=vTT 



vb 



vb 



(2tt) 3 T 3 

The fluid equations then become 

c 3 T 

v w c 2 (ijf + n bg ) + v w c 3 (5T' + 8T bg ) + — (v' + v bg ) + ^r^i + 5TT T1 



v w c 3 (fi f + /4 ) + v w d(5T' + 8T b[ 



bgJ 



^( v ' + v ' bg ) +tl r, 2 + 8TT T2 



.01614T (41) 



2T K 1 

V W C 2 / 2\/ 

m 



2T 



c 3 



C 4, 



v w c 4 T , 



g \r- ■ r bgj 1 g v."- 1 ~ bg) ■ g 



v' + vD+vTT, 







(42) 



where for top quarks one should use c/ and Tf and for W bosons one should use c b 
and T b . 
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We caution the reader not to interpret the T simply as rates for processes, without 
including the coefficients which appear in the derivative terms of the fluid equations. 
For instance, the rate at which a chemical potential for top quarks decays is not 
F^if ~ T/110, but roughly ~ r M i//c 2 / ~ T/9, corresponding to a typical lifetime 
for a top quark before it annihilates with an anti-top of 18/T (the factor of two is 
because each annihilation destroys two top type particles). Similarly, the time it 
takes for successive small angle collisions to randomize a particle's velocity is roughly 
c 4 /3r<, ~ 10/T for top quarks. 

Collisions between massive species and light species appear in the fluid equations 
of the light species with opposite sign. Since the light species are treated as being 
at a common temperature and velocity, no ST^Tt or v bg T v appears (see Appendix 
A). (If we treated the background species separately, terms order a 2 s \nl/a s would 
damp the difference between background species temperatures, velocities, and chem- 
ical potentials, efficiently forcing them to equal.) The background chemical potential 
is damped, but only by inelastic processes which enter at 0(o? s In l/a s ). The back- 
ground fluid equations are 



N w (n h Tn lh + ST b T Tlb ) 



Nw(^n2b + ST b r T2b ) 
N t v f TT vf + N w v b TT vb (43) 

which simply state that the particle number, energy, and momentum lost to the 
massive species via collisions are taken up by the light species. Here I]c 4 = 78c 4 / + 
19c46 = (87.25)47r 2 /30, which is the heat capacity of the light degrees of freedom, 
and N t = 12 and Nyy = 9 are the number of degrees of freedom of top quarks and 
W bosons. The collision rate T^g is the average over all background species of the 
particle number destruction rate. 

Let us estimate the importance of /i bg ] consider a very thick wall so that all 
derivatives may be dropped, and for simplicity ignore W bosons and ST and v. Then 
there are two fluid equations to consider, 

NhgT^gfihg = NtfitT^ T^fit = v w ^(m 2 )' 

The friction on the wall depends on the chemical potential of the top quark, which is 
+ Hbg- In this approximation, it is 

We see that \i bg is important when N t T^i > N bg Y ^ bg , or when the total rate at which 
background particles convert into top quarks exceeds the total rate at which they 



C ^W^ b g + C 3( V U,ST b ' g + l ^) + NkgT^gfihg 



3 

N t (v f r„ lf + 5T f r Tlf ) + 

Tvi, 



3 ' 
8TfT T2 f) 

~^^b 9 + -^-{^T bg + v w Tv bg ) 



18 



annihilate via inelastic processes. If these inelastic processes were very inefficient, the 
friction on the wall could be significantly enhanced. Unfortunately it is quite hard 
to compute as there are many diagrams, all with 5 particles. The diagrams 

involving several gluons interfere and are not separately gauge invariant. While the 
diagrams are at a high perturbative order, they gain infrared log enhancements when 
there is a t channel exchange and one incoming gluon is soft. The diagrams involving 
gluons have very large color factors because the gluons are in the adjoint represen- 
tation, and Bose statistics also make these diagrams large. We have made a crude 
estimate of these rates and find that Nb g Tb g is significantly larger than N t T^if, so we 
will neglect the background species chemical potential in what follows. (This conclu- 
sion sounds like an invalidation of the perturbative expansion in a s , but this is not 
so: N bg is several times larger than N t , and this makes up for the diagrams being at 
a higher perturbative order.) 

In this approximation the background fluid equations become rather simple 



c±{v w 5TL H = Nt^fT^f + 5T f T T2 f) + Nw^Y^b + 5T b Y T2b ) 



, Tvi 

bg 

^(8Tl g + v w Tv' bg ) = N t VfTT v f + N w v b TT vb , Hg = (44) 



Integrating these expressions, together with Eq. (22), reproduces Eq. (pip . The 
derivative term on the lefthand side is exactly what appears in the fluid equations for 
the perturbed species, which allows us to eliminate the background perturbations in 
the fluid equations by direct substitution. The final form of the fluid equations for 
the top quarks is 

C 3 f T 

v w C2fV' f + v w c 3f dT' f H J — v' f + HfTfAff + 5T f Y T i f f + [i b T^ lfb + 5T b T Tlfb = F lf 

C* A f I 

v w c 3f [i' f + v w c 4f dT' f + —^-v' f + fi/T/jzff + 5T f T T2ff + n b T ^fb + $T b T T2fb = F 2f 

°f»' f + C f5T' f + v ^L v ' f + VfT T vff + v b TT vfb = (45) 

where 





F 2f - 


v w c 2f , 2\l 

2 K) 


1 - 1 Ml/ + ^ 1 M2/ 
Zv c 4 


T Tlff 


r , N t C 3 f 
- 1 Tl/ + „ ± T2/ 


1 M 2// - I 1 + ^ M/x2/ 
2-*i c 4 


T T2ff 


- (1+ „ JlT2/ 
Z^ c 4 


r N w c 3fr 

Z^i c 4 


TtI/6 


N w c 3 f r 

— „ J- T26 


r N w c 4fr 

»W - sr n 1 M2b 




" EC4 ^ 


1 vff - (1 + „ Jl vf 

c 4 


r vfb = 


N w c 4f 
LC4 
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The W fluid equations look the same with the replacements &«->/, Nw <-> iVf . 
The collision terms have become (weakly) coupled between the heavy species, now 
indirectly through their influence on the background. 

We also extract the behavior of the background temperature for future use: 

r T i NtjufT^f + 5T f T T2f ) + NwiubT^ + ST b T T2b ) 

Ec 4 (l/3-^) 1 > 

We see that the resistance to the wall's movement from the heating of the plasma 
becomes important as the wall approaches the speed of sound. The divergent behavior 
at the speed of sound signifies the breakdown of our linearization of perturbations. 



7 Delta function response 



In order to gain some intuition for the fluid equations we study the response to a 
delta function source. Consider the fluid equations for the top quark, ignoring for the 
moment the change in the background so there is no coupling to the W boson. We 
can write the fluid equations in a matrix form, 



AS' + T5 = F 



(47) 



where 



5 = 



( A 4 
ST 



v 



( 




C3V W 


C3 

3 




C3V W 


C A V W 


C4 

3 


V 


C'i 


c 4 


C4V W 


3 


3 


3 





( F 1 \ 




F = 


F 2 

Vo J 




/r„i 






















Since the fluid equation is linear it suffices to study the solution when the source is a 
delta function times some column vector F. To solve this we need the homogeneous 
solution, 

S = ^OiiXi exp(-Ajx) x = z + v w t 



where ctj are coefficients, Aj are the solutions to 

Det[-AA + T] = 



(48) 



and Xi is the vector annihilated by — A; A + T. 

To solve the problem with a delta function source, we write a general solution 
to the homogeneous equations on each side of the source and apply boundary and 
matching conditions. Because S should go to at large distances, only the positive 
values of A can have nonzero coefficients in front of the source and only the negative 
values can have nonzero coefficients behind. We write 



EA,>o«iXiexp-Ax x>0 
EAi<o a iXi exp -Xx x<0 



(49) 
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and determine a>i from A5(0 + ) — A5(0 ) = F, the matching condition across the delta 
function. This gives 

J2^gn(X i )a i A Xt = F (50) 

i 

which is solved by expanding A~ Y F in the eigenvalues \. 

We see that the solution consists of several tails, some in front of the source and 
some behind it. These tails model the transport of the perturbation around the source 
due to particle flow. 

It is interesting to know how far the particles spread and how asymmetric the 
spreading is, so we will very briefly investigate the roots to Eq. (f48|). This equation 
is a polynomial in A. The coefficients are very messy, but much of what we want to 
know is in the coefficient for A 3 , 

CaV w A 2 w 2\ 

-^-(3-^)^204-03) 

and the coefficient for A , 

r^(r /t irT2 — r^iYi) 

Their ratio gives the product of the A's, and its sign tells us how many of the tails 
precede the source and how many follow it, as a function of v w . We immediately 
see that, for a subsonic wall, one tail precedes the source and two follow it; but the 
sign of the coefficient for A 3 changes at v w = v s = l/v3, and all three roots then 
trail the source. No diffusion occurs in front of a supersonic wall, at least within the 
fluid approximation. It is thus quite important to the study of baryogenesis to know 
whether the bubble wall is subsonic or supersonic. 

We will also comment that, in the special case that v w = and the particle decay 
rate is very much slower than the scattering rate (so IY2,r„ >> F^i, Tti, T^), the 
root equation for A simplifies; it is approximately 

c i 

— -q-iy 2 a + r^r^r^i = o 

with roots A — > 00 (a non-propagating disturbance) and A 2 = QI^T^/cg, which are 
decay tails. The length of the tail is C3/ (3yT\T^i). By finding the small v w limit of 
the coefficient for A 1 and comparing the tails to the result of the diffusion equation 
we find that the diffusion length is D = c\/(Qc2T v ) [0, which for our value of T v is 
D = 2.7 '/T for top quarks and D = 5.5/T for W bosons. 

From the collision rates presented in the last section, we find the length of the 
tails at v w — is 7.4/T for top quarks and 18/T for W bosons. Top quarks do not 
spread very far, but W bosons do. However, both lengths are larger than the diffusion 
lengths of the species, so the journey is a random walk and not a free flight. 
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8 Wall Velocity with a Wall Shape Ansatz 



In the fluid approximation the equation of motion of the Higgs field, Eq. (Q), becomes 

- (1 - vlW + V+(<f>, T) + ^f^icim + c 2f 5T f ) + "f^ icim + c 2b 5T b ) = 

(51) 

This, together with the fluid equations and the equation for the background tem- 
perature, constitute a well posed set of equations for the shape of the wall. They are 
velocity dependent; we expect them to have a solution at a discrete set of velocities 
(hopefully one). However, they are nonlinear, so their solution must be numerical. 
For starters it would be useful to see how much progress we can make analytically. We 



have done so in an earlier paper [p2j] ; here we will improve that analysis by including 
the background temperature. 

Let us restrict the shape of the wall to an Ansatz, 

0M) = ^(l + tanh^^) (52) 

where v w and L are parameters to be optimized. If Eq. (^l|) were derived from a free 
energy F we would know how to proceed: we would solve dFj dL = J Eq. (fn\)d<p/dL = 
and dF/dv w = J Eq.(fn\)d(fi/dv w = simultaneously. Since our equation of motion 
is dissipative, there is no free energy which generates it; however, these constraints 
still have the right physical content. Noting that d<p/dL = —{z + v w t)<f)'/L and 
d(p/dv w = t(f> f , we guess that a good pair of constraints should be 

J (Eq.51)0'dz = J (EqlOD-^'tfe = (53) 

Indeed, these have sensible physical interpretations. From Eq. (|2|), we see that the 
first of these constraints is that the total pressure on the wall should be zero; if the 
total pressure were nonzero then the wall would accelerate, changing v w . The second 
equation is an asymmetry in the total pressure between the front and back of the 
wall. If it were not zero there would be a net compressive or stretching force on the 
wall, changing L. 

Next we must deal with the variation of temperature across the wall. Write the 
temperature in the broken phase as T (which is T a in the language of section [|) and 
the temperature at a position z as T(z) = T + ST bg (z). We will solve for ST bg using 
Eq. ( fl6| ) and the boundary condition 5T bg (z — > oo) = and correct (|5"ID to linear 
order in 5T bg . The correction is 5T bg dV^ / dT ~ AD(pT5T bg . 



The integrals (|53D for the constant temperature, equilibrium part of (^T|) can be 
performed; they give 

p(f> + V^(4>))4>' = V T ((po) ~ V T (0) = -AV T (54) 
□0 + W)]f0' = ^^-^[AV T + E} (55) 




Ctifa-l)!* + wP + ^H 2 - 79 +oln^ (56) 
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Note that the □</> term acts to stretch the wall (increase L) while Vr acts to accelerate 
and compress the wall. The coefficient 2.79 in the last term is the only place where 
our choice for the function Pit enters our computation. 
Next we must evaluate the integrals 

ADT J STbgM'dz, ADT J Z>T hg ^$dz 

for the background temperature, 

/y 2 f y 2 % 

Ntfaffif + C2f5Tf)—<fx/> (z)dz , J N t (c lf /j, f + c 2f 6T f )—-(j)(f)'(z)dz 

for quarks, and similarly for W bosons; to do this we Fourier transform the integrals 
to 

i y 2 Tiii , , dk 



J Nt (c lf fi f (k) + c 2f 5T f (k))^m-k)- 



et cetera. One complication is that is weakly z dependent; we approximate it at 
its value where 4>4>' is maximum. Evaluating (jxfi' by a contour integration, we find 

W{k) = y (l - *A;L/2)— csch— — — = (57) 

We determine p, and ST from the fluid equations; they can be written in a matrix 
form by writing 5 = [/if, STf, Vf, fib, 5Tb, Vb] T ; the fluid equations become 

AijSj + TijSj = F441 5% = RiSi (58) 

The coefficient matrices A and T and the source vector F can be read off from Eq. 
(fHf), and the form of R can be read off from Eq . fl46|) . 

Fi has one term ~ c\b which is weakly position dependent; again we approximate 
it as its value where cfxfi' is maximum. The fluid equations are then easy to Fourier 
transform, becoming 

ikdi + (A-^ijTjA = {A- l ) tJ F 3 H) ! ik5T bg = RiSi (59) 

Now denote the eigenvalues of (A~ 1 T) as A,, and write the matrix whose columns are 
the eigenvectors of (A^T) as %i so that 

{A 1 T) i jXjk = Xik^k 

(with no sum on k). Also define aj = (x 1 )ij^j an d Si = (x~ 1 )ij(A~ 1 F)j, which are 
the deviation from equilibrium and the source in the basis of eigenvectors. Multiplying 
the fluid equations on the left by x" 1 , we find 

5i = 77T#) ikSf bg = RiXij T^rW{k) (60) 
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The equation for 5Tf, g only determines its value up to a constant of integration, 
which in Fourier space is a delta function; the coefficient of the delta function is 
determined by the boundary condition for 5Tb g . We find 



ST. 



bg 



RiXij 



SjTC 

ik{\j + ik) Xj 



S(k) #'(*;) 



(61) 



where S(k) is the Dirac delta function. 

The contributions to Eq. (|53|) in this notation are 



and 



dk 

dk 
2^L 



fiXij 



S; 



Xj + ik 



ADTRiXij 



S; 



ik(Xj + ik) Xj 



Sj7T 5(k) 



fiXij' 



S; 



+ ADTR iX ij 



Sj7r -S(k) 



V(k)<W(-k) (62) 



(f)(f)'(k)z(f)<f)'(-k) (63) 



Xj + ik ' WJ \ik(Xj+ik) Xj 

The vector / above, which gives the force on the field from 5, can be read off from 
Eq. (0). 

How does transport enter these integrals? The A, are the inverse lengths of expo- 
nential tails. The no transport limit is the limit ij^whichthe A, are large. In this limit 
the ik in the denominator is irrelevant. Since <jxft'(k)<j)<f)'(—k) is real, the ik reduces 
the value of Eq. (|62D ; transport, by spreading out the perturbation, has reduced the 
friction on the wall. 



Eq. (|62|)- (|63|) may be solved using the following integrals: 



dk 



2tt A + ik 
dk 1 



<f>(f)'(k)<f>(f>'(-k) 



16 V 4 ; v 2 ' 



XL 
~3 



(64) 



2ttL X + ik 



(j)(p'(k)z(j)(j)'(-k) 



dk 



2wik(X + ik) 
dk 1 



<M/(k)<M/(-k) 



16 \ 

+(1 
-1 



(XLf 



XL 



ALtt, 



XL 



7T , 



(M/(k)z<M/(-k) 



2ixLik{X + ik) 
Where we have defined the integrals 

x 2 csch 2 x 



A 



A 



4 

(Eq.64) 



6 



(Eq.p|) 



96 



h{a) 



dx 



-oo x 2 + a 2 
Evaluating these by contours gives 

iia 



h(a) 



-°° x 3 csch rrcothx 
x 2 + a 2 



dx 



(65) 
(66) 
(67) 

(68) 



h(a) 
h(a) 



n(2ira) 



sin 2 a 



na cos a 



sin 3 a 



^ ((utt) 2 - a 2 ) 2 

3a 2 + (mi) 2 



2„2 



n=l 



((mr) 2 — a 2 ) 3 



(69) 
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Numerically, these expressions suffer from cancelling divergences near a = mr. 
For large values of a the integrals possess useful asymptotic series 

At small values of a it is useful to rearrange the infinite series into Taylor expan- 
sions, 

h(a) = ^-2-4$>(-) C(2n + 1) 
sin a ^ \7r/ 

_ 1 . 7ra 2 cosa . , . ( a\ 2n /3a 3 > . „ . „ A 

/ 2 a = — -3 1 + E nn + 1 " — C(3 + 2n) + C(l + 2n) (71) 

sin a ^ V 71 J \ 7T / 

which have radius of convergence 7r. 

This completes the evaluation of Eq. (^). These two constraint equations deter- 
mine curves in the space of v w , L and their intersections are self-consistent solutions 
for the wall shape and velocity within our Ansdtze and approximations. These curves 
are illustrated in Figure |2]. 



9 Solving the Equations of Motion 



It is also possible to solve (plD and the fluid equations numerically for a general wall 
shape. Although they constitute a nonlinear system of equations, they are linear in 
5; we can therefore solve for 8 as (nonlocal) functions of (f) and reduce the system to 
a single integro-differential relation for <p. We begin with Eq. (|59| ), but including the 
position dependence of F, so we must Fourier transform Fcfxj)' as a unit. Now write 

St = { X - X )ij{A- x ) jk F^' (72) 

We quickly find 

~ 5i = x -x^k (73) 

which can be inverse transformed into a convolution, 



6i(z) = ^ J dyG^y - z)^ 1 A- 1 F) 3 {y)<P < P\y) 

Gj(x) = sign(Xj)9(x sign(Aj)) exp(— Xjx) (74) 

Here G is the Green's function for fluid perturbations. Now given a spatial configu- 
ration for <p, say on a lattice of points, we may integrate numerically to find 5 and 
integrate again to find Tb g . This gives us the full equation of motion fl5*T|). If the 
equation of motion had been derived from a free energy, then we know that changing 
in the direction dictated by the equation of motion reduces the free energy, and 
is guaranteed to approach a minimum if one exists. We may hope that evolving (ft 
in the direction dictated by the equation of motion will lead us towards a solution 
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of the equations of motion, but we must be a little careful because there is only a 
solution at select values of v w , and because the wall has a zero mode. The naive way 
to evolve v w towards its correct value is to move it in the direction dictated by the 
total pressure on the wall, / (Eq.^)<p'dz. We then subtract a quantity proportional to 
(f)' from the equation of motion when we correct 0, as otherwise part of the correction 
we implement will just be a shift in the wall position. This approach should work 
below the speed of sound, but above the speed of sound the pressure may decrease 
with velocity (as the temperature elevation reduces). 

We have implemented this procedure numerically, and find that the wall does 
approach a solution of the equation of motion for subsonic starting velocity. Since 
the lattice on which is defined need only be 1 dimensional (as the wall is planar) 
it is easy to make the stepsize small enough that the results are robust; our results 
agree to 0.2 % when we change stepsize from 1/T to 2/T. We find that the wall is 
somewhat deformed from a tanh, appearing steeper on the side facing the symmetric 
phase and shallower on the side facing the broken phase. Part of this is due to the 
two loop terms in the effective potential, but it also turns out to be true when we use 
just the 1 loop potential, which produces a symmetric wall at equilibrium. We can 
understand this deformation as follows; the source for excess particles goes as <f)<f)', 
which peaks above <fi = (po/2. The excess particles are on average swept further up 
the wall towards large because of the motion of the wall with respect to the plasma. 
Hence most of the frictive force on the wall occurs at large 0, and stretches out the 
upper part of the wall. 



10 Results and Conclusions 

The wall velocity computed from the tanh Ansatz and numerically are tabulated 
for several effective potential parameters in Table [l] and the wall shape they give 
is compared in Figure We see that the Ansatz returns the velocity accurately, 
although at large Ap it gives a wrong wall shape. The column <fro/T is the Higgs vev 
after the phase transition has completed, accounting for the heating of the plasma 
due to the release of latent heat in the transition. 

There is always a subsonic solution within the approximations that we have made, 
because as one approaches the speed of sound from below the temperature is elevated 
by more and more. In all cases in the table with Ap = 0, the only solution is subsonic. 
This is because the liberated latent heat raises the temperature of the plasma, which 
inhibits the motion of the wall. Neglecting the background temperature, as we did 
in p2| , produces a supersonic wall, so the effect is quite important. The result that 
the wall is subsonic is robust in the sense that, if we have underestimated all collision 
rates by a factor of two, we still find only a subsonic solution. Also note that the 
velocity is quite weakly dependent on the Higgs mass; in the range < m H < 90Gev, 
.38 < v w < .46. Recall that this result ignores the Higgs particle contribution to 
friction and may underestimate the contribution of infrared W bosons, which can 
only make the wall velocity lower. However, top quarks were typically responsible for 
about 60% of the friction and 65% of the liberated latent heat, so we do not anticipate 
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m H (T = 0) 


Ap 


Ansatz v w 


Ansatz L 


numerical v w 


d>o/T 


.0204 


o 





.351 


29 


.365 


.987 


023 


34 


n 


356 


28 


374 


907 


.03 


50 


o 


.365 


26 


.392 


.757 


.04 


70 


o 


.377 


25 


.412 


.635 


.05 


81 


o 


.390 


24 


.428 


.562 


.06 


91 


o 


.401 


23 


.441 


.516 


.03 


50 


.1 


481 riiTi 


18 


.496 


.998 


.04 


70 


.1 


.497, .97 


15, 14 


.513 


.861 


05 


81 

Kj X 


.1 


508 91 


13 11 


.524 


777 


03 


50 


.2 


499 run 

. t: %J - J ■ J. LXJX 


15 


510 


1.12 


04 


70 


.2 


KIM run 

. '.J -L t.) . J_ LXJX 


12 


.524 


971 


.05 


81 


.2 


522 98 


11, 7.7 


.533 


.878 


.06 


91 


.2 


528 96 


9.6, 6.2 


.539 


.814 


.04 


70 


.3 


519 run 


11 


.530 


1.05 


.05 


81 


.3 


.527, run 


9.4 


.537 


.952 


.06 


91 


.3 


.534, run 


8.5 


.543 


.882 


.07 


98 


.3 


.538, run 


7.8 


.546 


.831 



Table 1: Wall velocity and thickness and Higgs vev after transition at several effective 
potential parameters. The second entries in the Ansatz velocity and length columns 
for some entries are detonation solutions; in some cases the Ansatz predicts runaway. 



that infrared particles will change our results a great deal. 

Including the background temperature has also reduced the dramatic stretching 
of the wall found in [Q. This stretching arose because the force on the wall from 



Sf is predominantly far back on the wall. In the limit where decay rates are fast 
compared to the wall passage time, 5f ~ <fxft', which peaks well above <fr — <f>o/2- The 
force on the wall from Sf depends on 5f<jxf)', which peaks even more strongly on the 
upper part of the wall. Transport compounds this effect as particles sweep up the 
wall because of its movement. The result is that most of the force from Sf is far back 
on the wall, stretching it out. However, STb g is largest (for subsonic walls) in front 
of the wall in the symmetric phase, so it tends to exert more force on the front of 
the wall, which compresses it. This partially compensates the stretching we found in 
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In all cases in the table in which Ap ^ 0, the Ansatz technique tells us that there 
is a supersonic solution; this is always very relativistic and sometimes accelerates 
without limit. While this result is not robust (it may be incorrect if friction from 
infrared particles is important), it is possible that the wall does become supersonic in 
these cases. However, the fluid approximation is not valid in these conditions, because 
the front part of the wall becomes quite thin. A more careful analysis shows that 
the compression of the front of the wall greatly increases the friction and prevents 
ultrarelativistic motion. This is discussed in Appendix C. 

In conclusion, if there is a gauge boson condensate in the symmetric phase, and if 
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neither this condensate nor infrared bosons impart substantial friction, then the wall 
becomes supersonic. On the other hand, if there is no gauge boson condensate in 
the symmetric phase, the wall is definitely subsonic just due to friction from thermal 
particles and the release of latent heat. 
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12 Appendix A: Scattering Rates 

In this appendix we discuss the computation of collision integrals. The collision 
integral appearing in the Boltzmann equation is, to second order in a, 

c ^ j^S£^ ms ' ms ^ +k - p, ~ k ' )vm (75) 

vm = hh(i ± ± Si) - fMi ± A)(i ± h) 

where the sum is over all four leg diagrams, p refers to the incoming particle, k is 
the particle it hits, and p' and k' are the outgoing particles; the legs are labelled 
as 1,2,3,4 respectively. M. is the scattering amplitude for the process. The /'s 
are population factors; the positive expression represents a particle being removed 
from the state with momentum p by a collision and is weighted by the population of 
particles in that state and of the state it collides with. The negative term in V[fi] 
accounts for particles scattering into the state. The factors 1 ± / for the outgoing 
particles arise from particle statistics; the ± is a — for fermions (corresponding to 
Pauli blocking) and a + for bosons (for stimulated emission). Weldon has given an 
excellent derivation of this expression from the discontinuity on the real time axis 
of the self-energy of the propagator fl32| . Extending his technique to include gauge 



particles in a general covariant gauge proves to introduce difficulties. The structure 
of the collision integral arises from the poles of propagators in a self-energy diagram, 
and the propagator D^ v = (g^ u — £,k^k v /k 2 ) /k 2 has extra on shell divergences for 



£ 7^ 0. Such infrared problems have been discussed by Braaten and Pisarski |33[], who 
have shown that it is sufficient to work in the gauge £ = 0. 

Let us begin by analyzing the population factors. Without loss of generality we 
may write /, = l/(exp a% ± 1), and noting that (1 ± /,) = fi exp dj we find 

hf 2 {\ ± / 3 )(1 ± U) - (1 ± /i)(l ± h)hf* = (e a3+a4 - e a ^)hhhf A 
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The equilibrium value of ai is j(Ei — v- pi)/T with a common v and T for all species. 
In this case, C[f] = by energy and momentum conservation. This holds for higher 
order graphs as well. In the Boltzmann approximation the collision integral is local, 
so it also holds for a common spatially varying temperature and velocity. 
Now let Oj = "f(Ei — pi ■ Vbg)/Tb g — Si, with Si small; then 

exp(ai + a 2 ) = exp(J^ j(Ei - $ ■ v)/T bg ) exp(-<5i - S 2 ) 

i 

~ expC£l(E l -p t -v)/T bg )(l-5 1 -5 2 ) 
to linear order in S, and hence 

v[f i ]~C£s)f 1 f 2 (i±f 3 )(i±u) 

where the sum on delta is over incoming minus outgoing legs. 

Now let us compute the 0(a 2 ) collision rates for top quarks. All t-channel di- 
agrams which contribute at order a 2 are shown in the first column of Figure |I[ 
There are also s-channel processes. These diagrams interfere, but luckily the cross- 
amplitudes are all 0(m 2 /T 2 ) and can be neglected. Top quarks can also decay di- 
rectly, with a matrix element O(g^), but this process is time dilated and phase space 
suppressed, which reduces its importance by 0{ra 2 / ET). It therefore contributes at 
order a w a y and can be neglected. 

When we compute the cross-sections of these diagrams using free particle propa- 
gators, we find that the t-channel processes lead to logarithmic infrared divergences. 
At finite temperature these divergences are cut off by the interaction of the exchange 
particle with the plasma, so the t-channel diagrams contribute at order a 2 \nl/a s . 
The exact thermal propagators are quite complicated, but if we are willing to make 
an error in the leading constant (but get the coefficient of the log right) we can 
make a simple approximation which renders the computation more tractable. For 
an exchanged (top) quark, the correct dispersion relation is very similar to that of 
a massive particle with m 2 = g 2 T 2 /6 EM. This correction is of order the (space de- 



pendent) contribution from the physical mass, m 2 = y 2 <ft 2 /2, which we will therefore 
neglect, making the collision integral </> independent. For the bosons the dispersion 
relations are quite complex. For space-like momenta, such as occur for t-channel 
gluon or W boson exchange, the longitudinal components of the propagator are De- 
bye screened, and the transverse parts Landau damped, below the plasma mass, which 
is m 2 = 2g 2 T 2 for gluons and m^y = (5/3)g^T 2 for W bosons f35fl . To approximate 



this effect we modify the denominator of the propagator, making it t — m 2 , where m 
is the appropriate plasma mass of the exchange particle. This approximation is quite 
rough; it should yield the correct leading log coefficient, but not the correct leading 
constant. Hence we can only work to leading log accuracy. We therefore drop s chan- 
nel processes, which do not produce logarithms, set u = —s (which is legitimate in the 
leading log approximation), and keep only the st/t 2 dependence of annihilation and 
absorption-reemission processes and the s 2 /t 2 dependence of scattering processes. 

We will also systematically drop terms of order m/T, which allows us to treat 
the outgoing lines as approximately massless. In all cases the collision rates are 
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dominated by thermal particles, and the infrared log divergence arises primarily when 
the exchange momentum is between T and m, so this approximation is justified. 

Top quark annihilation rates 
Now we will compute a simple annihilation diagram, quarks go to gluons with 
matrix element (in the leading log approximation) ~ — (64/9)^st/(t — m 2 q ) 2 . The 
collision integral, integrated over d 3 p/ (2n) 3 T 2 , is 



I 



(2n) 3 T 3 2E p {2nf2E k 
d 3 p' d 3 k' -st 



2n) A 5\p + k-p' -k')V[fi] (76) 



(2ir) 3 2E p , (2ir) 3 2E k , (t - m 2 q )< 
In the leading log approximation the energy transfer is small, so we may approximate 

v = f p f k (i + / P 0(i + AO - f P fk(i + f P )(i + fk) 

where the 1 + / use the Bose population function and the / use the Fermi population 
function. In this approximation, V is independent of p' and k' and the integrals over 
p' and k' are Lorentz invariant and may be performed in the center of mass frame. 
The integral over k' performs the three spatial delta functions, so that k' = —p'. The 
integral over p', making a small mass approximation for t, is 

p' 2 dp'da p , (2p.k)2pp'(l-co S 9') 1 2p-k 

-27cd{2E p -2E p/ )—— — — ~ in— — (77) 



(2n) 3 2E p ,2E k , y p ' p ' {2pp'{\ - cos 9') + m 2 ) 2 8vr m 



to leading log accuracy. We now perform the remaining integrals in the plasma 
frame. The argument of the log contains the only dependence on p or fc; 2p ■ k = 



2\p\\k\(l — cos9) + 0(m 2 ), with 9 the plasma frame angle between p and k. The 



>r 2 

remaining integrals are 

64(^ 1 r p 2 dpdQ p k 2 dkdipd cos 9 



f p dpd\l v k dkdipd cos V . , ^ ^ . ^ /rr[ , , . . 



fMl + /p)(1 + ft) 'n 2pk{1 ~ 2 COS9> (78) 

m q 



The integral over (p z + k z )v vanishes when we integrate over dQ p , although (only) it 
will contribute when we integrate over p z d 3 p/T 3 . The integral over 9 gives approx- 
imately 2 \n(Apk/m 2 ), again dropping a constant and keeping only the leading log. 
The remaining angular integrals are all trivial. The energy integrals are dominated 
by thermal energies where we are justified in approximating p = E p , and we get 

f , ro , / i 7 \ cm //Til i ^pk 



18tt 5 T 3 
Using the approximation 



/ pk{2^ +(p + k)5T/T] In -^-/ p / fc (l + f P ){l + f k )dpdk (79) 



Jp n \n^f p (l + f p )dp~Hn + l/2) Jp n f p (l + f p )dp (80) 
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which is justified asymptotically and is reasonably accurate already at small values 
of n, we perform these integrals and get 



9tt 3 V 16 m 2 4 4 m\ J K ' 

The integrals over E p d 3 p and p z d 3 p are performed similarly and introduce no further 
complications. 

Scattering Processes 

Scattering diagrams are somewhat more complicated because the sum over chem- 
ical potentials and energies involves both incoming and outgoing particles. Con- 
sider the t channel gluon exchange diagrams. The amplitude squared and summed 
over outgoing states is ~ A0gjs 2 /t 2 for a top quark scattering off a gluon and = 
(5/6)32g 4 (s 2 + u 2 )/t 2 ~ (160/3)g 4 s 2 /t 2 for a top quark scattering off a quark (when 
it scatters from another top quark, 8 — as the other top is at the same tempera- 
ture and velocity). We have replaced u — > — s which is legitimate in the leading log 
approximation. The sum over perturbations 8 — 0/i + (E p — E p /)ST + (p z — p' z )v 
so these diagrams do not damp chemical potential. Also, the integration measure 
is symmetric under p <-> p' and k <-> k', but the integrand is antisymmetric, so the 
diagrams do not contribute to the first fluid equation. Also, the integral p z d 3 p gets 
no contribution from ST because the integrand is invariant under parity transforma- 
tion pi — > —pi, with pi = {p, k,p', k'} and the integral over E p d 3 p similarly gets no 
contribution from v. Hence we need only two integrals, 



A s 



2 



pkp'k' T 4 (t — m 2 ) 2 



(2nrS(p + k-p'-k')f p f k (l-f pl )(l±f k ,)l [ E ^ P _lf (82) 



where A = 40gf for the scattering off a gluon, and A = (160/3)g 4 for the scattering 
off a quark and of course the correct population densities for bosons and fermions 
ought to be chosen. We use the symmetry of the integration measure and V[f] under 
p <-> p' to rewrite E P (E P — E p >) — > (E p — E p /) 2 /2 and write it in a Lorentz invariant 
form as [u- (p — p')} 2 /2, where is the unit vector in the time direction of the plasma 
frame. Similarly, p z {p z — p' z ) — > —t/6 + [u • {p — p')] 2 /6- 

Writing the center of mass frame angle between u and p as [3, the integrals over 
p' and k' give, in the center of mass frame, 

A I (\u\\p\ sin/5) 2 In 

8n\ H(\u\\p\smP) 2 + 2p 2 }\n^ 1 } 

A quick calculation gives 

(-^p -^k) plasma 



I u\ cm \p\ cm cos (3 = u-p = u- (p- k)/2 
and 

plas 

|^| cm | P| cm. ~Z 
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so, in terms of the plasma frame angle 9 between p and k, 



is 



A j E p E k (l + cos 6)ln(2E p E k (l - cos 6) /ml 



16n 



\E p E k (3 - cos9) \n(2E p E k (l - cos9)/mf 



»4) 



Recall that A = AOg* (A = (160/3)gf) for scattering off a gluon (quark). We can now 
do the remaining integrations as above, using the same leading log approximation as 
in the previous section, the result is 



5 

7T 



6 mi 




20 2 2 25T 2 / 1/2 
S^ln^T 1/2 



i5) 



for scattering off a gluon and a quark, respectively. Bose statistics have made scatter- 
ing from a gluon more important than from a quark, even though the matrix element 
squared is smaller. 

Results for top quark diagrams 

We now tabulate the contributions of the collision integrals to the decay constants 
T for the top quarks. 

The contributions of two annihilation diagrams, both with matrix elements ~ 
-(64/9)st/(* - m 2 ) 2 , sum to 



32oj 
9tt 3 



(9C|/16)ln(9T 2 /m2) r> 

(21^/16) ln(15T 2 /m 2 ) I> 

(21C 2 C 3 /16)ln(15T 2 /m 2 ) T T1 

(135(2(4/64) ln(21T 2 /m 2 ) + (49(1/32) ln(25T 2 /m 2 ) T T2 

(45C 2 C4/64)ln(21T 2 /m 2 ) I\, 



56) 



The scattering diagrams when summed yield 

' C|ln(25T 2 /m 2 ) 
C 2 ln(25T 2 /m 2 ) 



25a; 
^3 



T 



Tt2 



Absorption and re-emission of a gluon, with matrix element 
m 2 ) 2 , contributes 



17) 



-(64/9)st/(t - 



at 



135C2C4 ln(21T 2 /m 2 ) - 98(1 ln(25T 2 /m 

36tt 3 " 1 90C 2 C4ln(21T 2 /m 2 ) -98C|/31n(25T 2 / 



2 ) 
q 

'm 2 



Tt 2 

r„ 



(88) 



These collision integrals were evaluated in analogous manner to the gluon exchange 
scatterings; after integrating over p' and k! one obtains (16g 2 /97r)(7?;pcos/3) 2 — > 
(16^ 2 /9tt)(1/4)[(p- k') -u] 2 for the contribution to 5T and (16g 2 /27n)[(-fvpcosp) 2 + 
p 2 } -> (165f 2 /27vr)(l/4){[(p - k') ■ u} 2 + 2p ■ k} for the contribution to v. 



Evaluating ( jSq) - 
3 2 T 2 /6, we find 
r w/ = 0.03765T p 



0.12, 



m r 



2g 2 T 



;|88[) numerically, using a s 
= .008993T, r^/ = -01752T = TT Tlf ,"T T2 f 
By comparing the contributions from various diagrams one 



2 and m 2 = 
0.07172T, and 



finds that scatterings dominate T v , but IY 2 arises mainly from annihilations. 
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We should comment that we have left out one potentially important diagram, weak 
flavor changing scattering, which converts left handed top quarks to bottom quarks. 
Because of the linear Coulomb singularity, cut off by Debye screening, this diagram 
contributes at order a w , and numerical evaluation shows that its contribution to the 
decay rate of left handed top quarks is comparable to that from the annihilation 
processes we have considered. However, it only affects left handed top quarks, and 
even if the rate were infinite it would just share their chemical potential equally with 
left handed bottom quarks, reducing the friction on the wall by 3/4. (The rate at 
which thermal top quarks rotate between right and left handed on the wall is smaller 
than the annihilation rate.) We will neglect this fairly minor effect here. 

W boson diagrams 

For the W bosons the dominant annihilation processes are t channel conversion 
to a gluon by a quark and W - gluon fusion to a quark- ant iquark pair. Summing 
over generations, flavors, colors, and particle-antiparticle, we find the matrix element 
is ~ —2^g 2 g 2 w stj{t — m 2 ) 2 for each process. (Again we only consider the leading log 
so s ~ — u.) We will also include order processes, which are double W fusion 
to fermions, with matrix element ~ — 18st/(t — m 2 ) 2 , W scattering from a fermion, 
with matrix element 120g^s 2 /(t — m^) 2 , and absorption re-emission, with matrix 
element ~ — 18g^st/(t — m 2 ) 2 . Scattering from another W boson does not contribute 
to the decay rates we consider because the sum of particle number, E, and p over 
all particles is zero. These processes do contribute to the damping of higher order 
perturbations, however, which may help to thermalize infrared W bosons. 

We will neglect annihilation to and scattering from the Higgs doublet because 
the matrix elements are more complicated, because the diagrams introduce infrared 
problems, and finally and most importantly because there is only one doublet, so its 
contribution is much smaller than the 12 fermion doublets. 

The collision integrals are completely analogous to those discussed above, so we 
will only present the results. For semi-strong annihilation, we find 



6 



7T U 



Ct.MnT < 



(9£/16) ln(9T 2 /m 2 ) 

(21C 2 C 3 /16)ln(15T 2 /m2 

(21C 2 C 3 /16)ln(15T 2 /m2 



(89) 



(135C 2 C 4 /32) ln(21T 2 /m 2 ) 
(45C 2 C4/32)ln(21T 2 /m2) 

For doubly weak annihilation, we find 

(9Cf/16)ln(9T 2 /(m 2 » 
(2IC2C3/I6) ln(15T 2 /(m 2 » 
(2IC2C3/I6) ln(15T 2 /(m 2 » 

(135C 2 C4/64)ln(21T 2 /(m 2 )) + (49(1/32) ln(25T 2 /(m 2 )) 
(45C 2 C 4 /64) ln(21T 2 /(m 2 » 

where ln((m 2 )) = 3 ln(m 2 )/4 + ln(m 2 )/4. 

For collisions from quarks and leptons, we find 



9 
2^3 



a 2 T 



Tt2 

r„ 



(90) 



15 
2^3 



<X„T 



2C 2 2 ln(25T 2 /m 2 w ) 
2C 2 2 ln(25T 2 /m 2 v ) 



Tt2 

r„ 



(91) 



33 



and just as in the case of quarks, scattering diagrams contribute equally to the IY2 
and r„. For absorption and re-emission from fermions, we find 

9 2 J (135C2C3/32) ln(21T 2 /(m 2 » - (49(1/16) ln(25T 2 /(m 2 » T T2 
4tt 3 w \ (45C 2 Cs/16) ln(21T 2 /(m 2 » - (49(1/48) ln(25T 2 /(m 2 » Y v ^ 

The thermal mass of a left-handed lepton is mf = 3^T 2 /32 plus a small hy- 
percharge correction which we neglect. The results, using a w = 1/30, are r Ml = 
0.00521T, r M2 = 0.01012T = T T1 , T T2 = 0.03686T, and r„ = 0.01614T. Annihila- 
tions dominate even r,,. 



13 Appendix B: Critique of the Fluid Approxima- 
tion 

Is the fluid approximation any good? 

The discussion can be broken into two parts: does the fluid approximation model 
properly the energy dependence of Sf; and does it oversimplify the anisotropy of 8fl 

We begin with the energy dependence, using the tools developed in Appendix A. 
Note first that top quarks decay fairly quickly (we have seen that the tails around a 
source have a length of about 5/T), so a particle excess decays before it is transported 
off the wall. Hence the friction from one excess particle is roughly its force on the 
Higgs field at the point where it was created times its lifetime, oc r/E. We must 
determine the behavior of r. 

The annihilation rate of one particle can be computed by putting a bump in the 
population function / at a specific energy and examining the collision integral. The 
contribution to 5 due to the bump is — 1/ f' Q times the bump, and f' Q = — / (1 — fo) for 
a fermion, so the population factor in the collision integral is /2(l±/3)(l±/4)/(l — fi) 
times the bump. The full energy dependence of the collision integral is this expression, 
the energy conserving delta function, and the 1/E p prefactor. The argument is that 
1 — fi is quite weakly energy dependent (going from 1/2 at E — to 1 at large E 
where most of the particles are) and for thermal outgoing particles the integral over 
outgoing states depends on E only through ln(2p • /c/m 2 ), a weak dependence. The 
dominant dependence is then the 1/E prefactor, which cancels the \jE strength of 
the particle's influence on the Higgs field. The argument is strongest for thermal 
particles, but as they dominate the total friction from top quarks, this should suffice. 

Let us check the performance of the fluid approximation for top quarks by ne- 
glecting spatial derivatives (ie assuming the particles really do not leave the wall) 
and comparing the friction from the fluid equations to the friction we would find if 
we had only included a fi term. Including only a ji term, we find 

liY^i = v w c\mm! friction = ci/imm' 

The fluid approximation gives 



H 

ST 



v w mm' 



C\Yt2 — c 2 Tti 
C2IV1 - Cil^ 2 



friction = (ci/i + c 2 ST^)mrri 
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Plugging in the values for the T's, the friction is v w m 2 m' 2 times .137 in the first case 
and .143 in the second case; using the full fluid approximation has only changed our 
estimate of the friction by 4%. (If the relaxation time approximation were valid the 
value would have been different by 27%, so the argument must have some validity.) 
Since we have difficulty trusting our evaluation of the collision integrals to better than 
50% (mainly because of the leading log approximation), and since the fluid equations 
we use should account reasonably well for transport, there seems little reason to 
improve the fluid approximation for top quarks, except by improving the precision of 
our evaluation of the collision integrals. 

The argument for the fluid approximation also works for thermal W bosons, where 
1 + / ~ 1; but for soft W bosons, the population term in the collision integral 
is / 2 (1 ± / 3 )(1 ± / 4 )/(l + /i) ~ / 2 (1 ± / 3 )(1 ± U)E/T, which cancels the l/E in 
front of the fluid equations; apparently, the decay rate does not rise as we lower the 
energy. Evaluating the annihilation and scattering rates of the infrared W bosons 
is further confused by the appearance of loop diagrams (hard thermal loops) which 
enter at the same parametric order as tree level effects. Also we should include W 
boson scattering from other W bosons, which does not contribute to the decay rates 
we have computed but does contribute to the rate at which infrared particles are 
thermalized. We will not attempt to treat this problem here, but will only remark 
that the fluid approximation appears to be a very naive treatment of the W bosons. 

Next, we will discuss whether it is acceptable to treat the anisotropy of the dis- 
tribution keeping only the lowest term, v -p — pViYu(p), rather than including higher 
angular moments p l Yi m (p), where Y tm is a spherical harmonic. The source of 5f is 
isotropic, so v only arises out of spatial gradients of 5f. The decay of v is character- 
ized by the diffusion length D ~ 2.7 /T for quarks (and 5.5/T for W bosons), which is 
generally smaller than the thickness of the wall. Relative to the other perturbations 
its amplitude is down by D/L. We expect perturbations at high angular moments 
to be down by (D/L) 1 ; as long as D is sufficiently less than L, we expect them to be 
less important than v. To confirm this we can extend our treatment beyond the fluid 
approximation to account for rank two tensor deviations. Setting T = 1 for simplicity 
(we can restore it by dimension counting) and writing 

f- 1 = exp(E-5)±l 

5 = fi + EST + piVi + E 2 e + Ep^ + (piPj - 5 ij p 2 /3)e ij 

we find extended fluid equations; 

c 2 fi + c 3 5T + c 4 e + ^diVi + y^e* - ciinm = -T^fx - T T1 5T - T el e (93) 
c 3 /i + c 4 ST + c 5 e + ^diVi + -^di€i - c 2 mm = -T^H - T T 2^T - T e2 e (94) 
c 4 /i + c 5 5T + c 6 e + — + y^e; - c 3 mm = -T^fi - T T3 5T - T t3 e (95) 

= -T vl v - r ell e, (96) 
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F v2 v - T ei2 ei (97) 
2(5 

(diVj + djVi - ^d k v k ) + ^(d^j + dfr ^-d k e k ) + 



J^K^ij + e 3i Y^ kk > = _r «J' e ^ ( 98 ) 

The matrix of T's is block diagonal in angular moments and symmetric. It contains 
10 distinct nonzero terms; its evaluation is the main impediment to including higher 
tensor moments. 

Only the isotropic moments \i, ST, and e are directly sourced, since only their 
time derivatives appear in the equations with the source terms cmm. The vector 
moments i>j and q are sourced by gradients in the isotropic moments. We expect 
their amplitudes to be ~ fi'/T ~ fiD/L. Similarly, the traceless tensor moment is 
sourced by gradients in the vector moments, d(Vy Although we have not calculated 
T ei j, we anticipate that it will be at least as large as the damping rates for the vector 
moments, because a particle's contribution to is erased by deflecting the particle 
by 45°, while it must be deflected by 90° to erase its contribution to i>j. Hence, 
should be ~ V{D/L ~ fiD 2 /L 2 . We conclude that as long as L » D, the neglect of 
high order angular moments should be justified. (In the opposite limit, D » L, we 
know that, for quarks, the high order moments are very important and lead to a log 
enhancement of the friction; so L » D is not only a sufficient condition to neglect 
high angular moments, but a necessary one.) In our case, D ~ 2.7 /T for top quarks, 
and we are in good shape. (Even for W bosons, where D ~ 5.5/T, we are in fairly 
good shape.) 

We wish to note that this conclusion relies on the fact that we are interested in 
the motion of one fluid against another, and not in the dissipation of the motion 
of a fluid as a whole, or of a fluid of one species. Let us explore this case briefly, 
to understand the physics of the extended fluid equations. When tracking a one 
component fluid, or the average over species of a several component fluid, energy and 
momentum conservation ensure that r„i, IY2, IVi = T^ 2 , T«2 = r ei i, and IY3 = r e2 
vanish identically, so ST and V; L are not directly dissipated. If in addition there is a 
nonzero conserved charge density for which /i is the chemical potential, then Y ^ and 
r e i = T M 3 also vanish, and /1 is not directly dissipated. It is still true that ~ ViD/L, 
but now its presence is important, as it is the main source of dissipation in the system. 
Indeed, if the perturbations are slowly varying, then we can neglect compared to 
d^j, et cetera, and solve for the high order perturbations in terms of the lower order 
ones. We find 

€ij = l^~^ iVj + 9jVi ~ ~^ LdkVk ^ 
'c 5 Vi + c 5 diST + c 4 dijj) = 777^— (c 4 - ^-)difi 



£ = —(^OiVi + C b ST + C 4 /i) 
1 e3 <J 
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all plus 0(5"/T 2 ). Here we have used Eq. fl96|) in the equation for e^. We can also 
use Eqs. (f?3|) and (^) to simplify the relation for e, but here we should note that 
the coefficients on diVi and 5T only equal in the massless limit. If the mass is zero, 
then e = 0; otherwise it is — 0(m 2 )diVi/T e3 . 

Substituting these three quantities into the fluid equations for the undamped 
species, we find 

c 2 li + c 3 6T + —diVi = c x mm H — d fx H 

c 3 fi + c 4 5T + ^-diVi = c 2 mrh + ° 4 ^ 5 C3 ° 5 d 2 /i + - <9^ 

O yC4l t i2 1 e 3 

c 3 c 4 c 4 . c\c^-c 3 c\ 0{m 2 ) 
—d i Li+—d i 5T + —v i = — — <9i/i + — d^ 

O O O yC4l e j2 J- e3 

2c 2 25- • 

+ 225T~^ ( ^^ + 9jVi ~ ~Y dkVk) (") 

which are the linearized relativistic Navier-Stokes equations. The perturbations e, e*, 
and 6jj have caused bulk viscosity, thermal resistivity, and shear viscosity, respectively. 
The bulk viscosity vanishes in the ultrarelativistic (m << T) limit; the thermal 
resistivity vanishes in the absence of a nonzero conserved charge density. (In this 
case heat flow is only resisted by higher derivative terms.) 



14 Appendix C: Runaway Wall 

According to the tanh Ansatz, the wall is capable of becoming ultrarelativistic without 
contracting to a small plasma frame thickness. We will explore how this arises and 
why it is unphysical. 

First consider the equation of motion in the tanh Ansatz, Eq. (|5^) and particularly 
Eq. (|55D . We see that the derivative term 0<p stretches the wall, but that it becomes 
ineffective at large velocity. One would expect, then, that the wall becomes extremely 
compressed. But the frictive term Eq. fl63|) stretches the wall. For a thin wall, where 
the particles have little time to decay, and at large velocity, where they simply sweep 
backwards up the wall, we find \x oc 2 and the ratio of stretching to friction is 
11/12. When H < 5AVr/6, then if there is enough friction to stop the wall from 
accelerating then there is enough to prevent it from contracting to the regime where 
we can neglect particle decay, so the wall cannot become thinner than about L ~ 6/T. 
When A F ^ 0, the system tends to supercool quite heavily and S tends to be small 
compared to AVr] the wall is then susceptible to runaway, according to the tanh 
Ansatz. 

Note that most of the frictive force on the wall comes on the upper part near the 
symmetric phase. The friction tends to make this part of the wall very thick, as we 
see explicitly when we solve for the wall shape. What the tanh Ansatz does is force 
the front part of the wall to be as thick as the back, which is unphysical. 

Let us attempt an analysis without a wall shape Ansatz. Consider the wall propa- 
gating at a steady, very relativistic speed, say jv w > 10. To determine how abrupt the 
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wall is, we will integrate the equation of motion times ', starting in the symmetric 
phase and going up to the point = <f>i\ if the wall is in a steady state, 

/_^ = Vt(0O*)) + friction)0' dz = (1 - vl) J <f>'cf>"dz = (1 - t£ (100) 

Because 1 — is very small, the wall is rising very rapidly where = 0i unless 
Vt{4>i) + / friction^' dz is almost zero. Let us examine whether this condition is ever 
satisfied. 

First note that the friction is never negative anywhere for a monotonic wall. Now 
Vt is positive at small 0, as otherwise the phase transition would have already pro- 
ceeded by spinodal decomposition. So there will be a section at the front of the wall 
where 0' is large, ie a section which is very abrupt. To determine where this section 
ends, we need to find J friction^' dz on this section of wall. 

This is easy in the fluid approximation. As the wall is thin, we can neglect decays, 
and the fluid equations can be integrated. The friction is, in the notation of section [8], 
fi(A~ 1 )ijFj(pf/8. Adding this to the effective potential, we get a curve whose second 
minimum lies above for the effective potential parameters in Table [j] with Ap = 0.1 
or 0.2, but whose second minimum is below for Ap = 0.3. In the former case, the 
whole wall is abrupt and feels a net backwards force; the wall will then slow down 
until it is not abrupt and no runaway occurs. In the latter section of the wall 

feels a net forward pressure, and accelerates without limit. According to the fluid 
approximation, the wall will run away in these cases. 

When the wall becomes this thin, the fluid approximation is definitely inaccurate- 
as we have argued, it only makes sense when the wall is thicker than the diffusion 
length. The free particle approximation should be appropriate, however, and for an 
ultrarelativistic wall we can make an accurate expansion in large 7. 

The total pressure on this section of wall from one species is, in the free particle, 
1 loop approximation (|5|, [0|) 

dpz , _ fYZ — 2^ f d2 P^ P z \ , 

2vr [Pz V Pz m) J (2tt) 2 y^-Texp^G, - vp z )) ± 1 + 

dPz 2 f d 2 p± p z 1 

2tt Pz J (2tt) 2 ^2^exp(/3 7 (p - vp z )) ±1 1 ] 

[°° dp z , — r d 2 P± Pz 1 

Jo 2vr l P* + VP* +Tn n (27r) 2 ^2 + ^ + m 2exp(/5 7 (v^+^ + ^))±l 

where m is the particle mass when = 0i. 

Taking 7 >> T/m so that exp(— jm/T) << m 2 /T 2 , we find the friction is 

2 T 2 m 3 Tln2 m 2 T 2 

+ - I5V + °(™ 4 K7r/m)/7 2 ) (102) 

from a fermion degree of freedom and 

2 T 2 s T / 5 ln(27T/m)\ m 2 T 2 n , . , . w 9x . . 

+ .1347 + y — 1 ' + O m 4 n( 7 T/m) h ) 103 

24 4vr 2 7 I 18 3 J 96-f 2 V U 1 ' J K ' 
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from a boson[]. Subtracting from these the corresponding 1 loop contributions to the 
effective potential, 

m 2 T 2 m 4 „ , , s . m 2 T 2 m 3 T , . 

~ (In vrT/m - 7^ + .75 , 104 

48 32vr 2 ^ 7 ' ' ' ' 24 12vr V ; 

gives the friction from particles on the abrupt part of the wall, at 1 loop. This 
friction is much larger than the prediction of the fluid approximation. Adding this 
expression to Vt produces a function which is positive for all 0i > 0, for all values of 
the effective potential parameters we have considered; in fact it comes quite far from 
having a second zero. This means that there is no value of where the abrupt part 
of the wall can end; it always has a net backwards force on it, even when (f>i = <f)Q. 
It is simply impossible to have an ultrarelativistic, steady state solution to the wall 
shape. As the wall becomes very fast, the front of the wall always compresses and 
sustains enough friction to prevent further acceleration. We cannot determine the 
velocity of the wall, but we can say that it must be slow enough that the free particle 
approximation is inaccurate, which ensures that ^v w cannot be much greater than 1. 
Note that it was not necessary for the condensate responsible for Ap to exert any 
friction for this to be true. 
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Figure 1: All annihilation, scattering, and absorption re-emission diagrams considered 
in the text 



Figure 2: Plot of the solutions of Eq. fl53l) for the case At = 0.03, A F = 0. The solid 
line is the velocity constraint and the dotted line is the thickness constraint. 



Figure 3: Wall shape; at left, for At = 0.04, A F = 0.1, and at right, for At = 0.03, 
A F = 0. The solid line is the numerical value, the dotted line is the tanh Ansatz 
shape. The scale is the same for both; the full width of the pictures are 96/T and 
64/T. 
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